Spin fluctuations and unconventional pairing in KFe2As2 
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We study the relation between the spin fluctuation and superconductivity in a heavily hole doped 
end material KFe2As2. We construct a five orbital model by approximately unfolding the BriUouin 
zone of the three-dimensional ten-orbital model obtained from first-principles calculation. By ap- 
plying the random-phase-approximation, we obtain the spin susceptibility and solve the linearized 
Eliashberg equation. The incommensurate spin fluctuation observed experimentally is understood as 
originating from interband interactions, where the multiorbital nature of the band structure results 
in an electron-hole asymmetry of the incommensurability in the whole iron-based superconductor 
family. As for superconductivity, s-wave and d-wave pairings are found to be in close competition, 
where the sign change in the gap function in the former is driven by the incommensurate spin fluctu- 
ations. We raise several possible explanations for the nodes in the superconducting gap of KFe2As2 
observed experimentally. 



I. INTRODUCTION 

A number of interesting features of the iron-based su- 
perconductors have attracted much attention. Among 
them is the material dependence of the superconductiv- 
ity. Not only the transition temperature Tc, but also 
the form of the superconducting gap seems to be ma- 
terial dependent. For instance, although a number of 
experiments suggest that the gap is fully open on the 
Fermi surface on some of the arsenides, there are sev- 
eral experiments suggesting the presence of nodes in the 
superconducting gap LaFePOi"— , BaFe2(Asi_a:Pa;)2^"— , 
or in LiFeP— . Also in the electron-doped 122 systems, 
Ba(Fei_a;Ma;)2As2(M=Co,Ni), the presence of gap nodes 
has been suggested^i^ 

For LaFePO, theoretical studies suggest the presence 
of nodes on the electron Fermi surface around the wave 
vector (tt, 0)/(0, tt) in the unfolded BriUouin zonejiSrJl 
From this viewpoint, the heavily hole doped end mate- 
rial KFe2As2 is of special interest because in this ma- 
terial, the electron pockets present in most of the iron- 
based superconductors are known to be absent due to 
the low position of the Fermi energyJ^ii^ Nonetheless, 
the material becomes superconducting, and interestingly, 
several experiments suggest the presence of nodes in the 
gapJi"— Since there are no appreciable electron pockets 
in KFe2As2, the nodes are likely to be on the hole Fermi 
surface. Theoretically, a possibility of d-wave pairing has 
been proposed recentlyi^i'^^ 

Another interesting observation in KFe2As2 is the spin 
fluctuations found in neutron-scattering experiments. In 
Ref. i23l . incommensurate spin fluctuation has been ob- 
served, and we have theoretically provided its explana- 
tion in terms of the multiorbital nature of the Fermi sur- 
face. 

Given this background, the aim of the present study 
is to understand the relation between the spin fluctua- 
tions and the superconducting gap of KFe2As2. We first 
construct a ten-orbital model for KFe2As2 and also ob- 



tain a five-orbital model by approximately unfolding the 
BriUouin zone. We adopt three-dimensional models in 
order to take into account the possible variation of the 
superconducting gap along the kz direction. We apply 
random phase approximation to these models and com- 
pare the calculation results at high temperature, which 
gives reasonable agreement, although the BriUouin zone 
unfolding is not exact. Then, we go down to low tem- 
perature using the five-orbital model. We first summa- 
rize our understanding on the origin of the incommen- 
surate spin fluctuations^'^ and also analyze its strength 
when the material composition is varied from BaFe2As2 
to KFe2As2. We finally discuss the superconducting gap 
when Cooper pairing is mediated by spin fluctuations. 
An analysis based on the linearized Eliashberg equation 
shows that s-wave and d-wave pairings are found to be in 
close competition, and the former is found to have essen- 
tially the s±-wave form in the sense that the gap function 
around the wave vector (0, 0) (where the hole Fermi sur- 
faces exist) and around (tt, 0)/(0, tt) (where the electron 
Fermi surface would have existed in a less doped mate- 
rial) has the opposite sign. This sign change is driven by 
the incommensurate spin fluctuation, which is due to the 
interband interaction between states around (0,0) and 
(tt, 0)/(0, tt). Since the electron Fermi pockets are ab- 
sent, the latter states are somewhat away from the Fermi 
level, so that this interband sign change of the gap func- 
tion should not be detected experimentally. As for the 
possible explanation of the nodes of the gap observed ex- 
perimentally, we raise several possibilities concerning the 
hole Fermi surface. 



II. BAND STRUCTURE AND FERMI SURFACE 

First we perform a first-principles band calculation for 
KFe2As2 adopting the experimentally determined lat- 
tice structure^! and using the QUANTUM ESPRESSO 
package;^ Then we exploit maximally localized Wannier 



2 



orbitals^S to obtain the ten-orbital model, whose band 
structure and the Fermi surface are shown in Fig. [T] (a) . 
There are ten 3d bands because there are two iron atoms 
per unit cell in the body-centered tetragonal (BCT) lat- 
tice structure. Around the X point, the electron Fermi 
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FIG. 1: (color online) Band structure and Fermi surfaces of 
(a) the ten-orbital model in the original Brillouin zone of 
the BCT lattice structure and (b) the five-orbital model of 
KFe2As2 presented in the unfolded Brillouin zone of the ST 
lattice structure. The dots in the left are the original first- 
principles band calculation, while the solid lines are the tight- 
binding model band dispersion obtained by using maximally 
localized Wannier orbitals. The horizontal cuts of the Fermi 
surface at fcz = and — tt are presented in the bottom 
panels of (a) and (b). Translating the left (right) panel of 
the five-orbital Fermi surface (b) by (tt, tt, tt) and superpos- 
ing it to the right (left) corresponds to the ten-orbital Fermi 
surface shown in the right (left), (c) The BCT and the ST 
lattice structure. The former (latter) unit cell with two (one) 
iron(s) is adopted for the ten (five)-orbital model. The figure 
has been produced using VESTA2I. (d) The Brillouin zone 
correspondence between the ten-orbital (black line) and five- 
orbital (blue line) models. The Z point of the second Brillouin 
zone of the BCT structure coincides with the M point of the 
unfolded Brillouin zone. 



surfaces, which are present in iron pnictides with fewer 
holes, are absent. Nonetheless, the band that gives rise 
to the electron Fermi surface for smaller number of holes 
sits very close to the Fermi level. Three cylindrical hole 
Fermi surfaces are present around the F point, which we 
name ai, 7, and 02 from the inner one to the outer. Also 
around the Z point there is a tiny hole Fermi pocket. 

For 122 systems, the Brillouin zone unfolding 
procedure^^ that adopts the reduced unit cell (with one 
iron) cannot be done strictly,® namely, the hopping inte- 
grals that go out from and come into Fel and Fe2 in Fig. 
[Ijc) are somewhat different. Nonetheless, if we adopt a 
reduced unit cell by neglecting the difference, we end up 
with an approximate five-orbital model.— We show the 
band structure and the Fermi surface of this model in 
Fig. [1] (b) in the "unfolded" Brillouin zone of the simple 
tetragonal (ST) lattice structure. The correspondence 
between the unfolded ST and the folded BCT Brillouin 
zone is shown in Fig. (Ud). As seen from this corre- 
spondence, the band structure of the five-orbital model 
can be refolded into the original BCT Brillouin zone by 
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FIG. 2: (Color online) (Top) The band structure of the five- 
orbital model of KFe2As2 refolded into the original Brillouin 
zone (solid red), and the band structure of the ten-orbital 
model (dashed black). (Bottom panels) Horizontal cuts of 
the Fermi surface for the ten-orbital (upper) and five-orbital 
(lower) models. The thickness represents the z^, xz/yz and 
xy orbital characters. 
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translating the bands in the second BCT Brillouin zone 
by a wave vector (7r,7r,7r). The comparison of the band 
structure between the original ten-orbital model and the 
five-orbital model in the refolded Brillouin zone is shown 
in Fig. [2] They coincide with each other in most of the 
portions, but there is some discrepancy along F-Z. 

In the middle panels of Fig. [51 we show the vertical 
cuts of the Fermi surface, where the thickness represents 
the strength of the z^, xz/yz, or xy — {X"^ — Y^) orbital 
character. All of the hole Fermi surfaces are very cylindri- 
cal, so that the two-dimensionality is strong. However, if 
we look at the orbital character, a2 and 7 Fermi surfaces 
have strong xz/yz and xy character, respectively, regard- 
less of fcz, while the ai Fermi surface is mainly composed 
of orbital around the Z point, while the xz/yz charac- 
ter becomes stronger around the F point. In this sense, 
ai can be considered as less two-dimensional compared 
to the other two-hole Fermi surfaces. 

In the bottom panels of Fig. [21 we show similar ver- 
tical cuts of the Fermi surface for the five-orbital model, 
refolded into the original Brillouin zone. The result in 
most portions is close to that obtained for the ten-orbital 
model, but the ai hole Fermi surface, which is a cylin- 
der in the latter, splits into a three-dimensional pocket 
around the Z point and a strongly warped cylinder. This 
split occurs between portions of the Fermi surface having 
strong and xz/yz character. Although this may seem 
like a large difference, the effect of this splitting is not so 
large as far as the quantities we focus on are concerned, 
as we shall see later. 



III. MANY-BODY HAMILTONIAN AND 
RANDOM-PHASE APPROXIMATION 



For the many-body part of the Hamiltonian, we con- 
sider the intraorbital U , the interorbital [/', the Hund's 
coupling J, and the pair hopping interaction J'. We 
consider orbital-dependent interactions;^ and the inter- 
action part of the Hamiltonian is given as 



The spin and the charge (orbital) susceptibility matrices 
are obtained as 



Xs{q) 



Xc{q) = 



i-5x"(g) 



1 + Cx0(<7)' 



(3) 



(4) 



where S and C are interaction vertex matrices. We re- 
fer to the largest eigenvalue of the matrices Sx^{q) and 
Xs{q) for iuJn = as the Stoner factor and the spin 
susceptibility, respectively. The magnetic instability is 
signaled by the Stoner factor exceeding unity. In the 
RPA calculation (where the self-energy correction is ne- 
glected), adopting the interaction values evaluated by 
first-principles calculation'^*' easily results in a magnetic 
instability even at high temperature, so we multiply all 
the electron-electron interaction by a constant factor /— . 

To analyze superconductivity, the Green's function 
and the effective singlet pairing interaction, 

V'{q) = lsx%q)S - lcx%q)C + ^{8 + C), (5) 
are plugged into the linearized Eliashberg equation. 



X(j)uj^{k - q)Gi,,i^{q - k) 



(6) 



The eigenvalue A increases upon lowering the tempera- 
ture and reaches unity at Tc. Instead of going down to 
Tc, which will be a tedious calculation, we perform the 
calculation at a fixed temperature and obtain the eigen- 
value and the eigenvector for s-wave and d-wave pairing 
symmetries. The comparison of the eigenvalue at a fixed 
(low) temperature for different pairing symmetries tells 
us which pairing state is actually realized in the system. 
The eigenfunction of the Eliashberg equation is referred 
to as the superconducting gap function. 
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IV. COMPARISON BETWEEN TEN- 
FIVE-ORBITAL MODELS 



AND 



JtivSifj, ■ Sin + ^ >^/ii/Cip-|-cj^^Ci,y-f-Cii,^ I , (1) 

where i denotes the sites and /i, v stand for the or- 
bitals. We apply random-phase approximation (RPA) 
to this model and obtain the spin and charge suscepti- 
bility matrices. Namely, using the bare Green's function 
Gim{k) {k = {k,iujn)) in the orbital representation, the 
irreducible susceptibility matrix is given as 

X?,.b./3.u('?) = EGw3(fc + '?)G^i.(fc) (2) 



We first compare the five- and ten-orbital models at a 
high temperature of T = 0.07eV. We take 16 x 16 x 16 
/c-point meshes and 128 Matsubara frequencies. The in- 
teraction reducing factor is / = 0.53. The Stoner factor 
calculated for this parameter set is as = 0.94 and 0.95 for 
the ten- and fivc-orbital models, respectively, in reason- 
able agreement. In Fig. [31 we show the spin susceptibility 
along the wave vector (0,0) to (tt, 0) in the unfolded Bril- 
louin zone. The slight difference in the Stoner factor re- 
sults in a difference in the peak value of the spin suscepti- 
bility, but the qx position at which the spin susceptibility 
is maximized is the same between the two models. It can 
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FIG. 3: (Color online) Comparison of the spin susceptibility 
between the ten-orbital (top) and five-orbital (bottom) mod- 
els along the wave vector (0,0) to (7r,0). T = 0.07eV with 
16 X 16 X 16 fc-point meshes. 



also be seen from this figure that the three-dimensionality 
of the spin fluctuations is very small. 

In Fig|4j we compare the superconducting gap for the 
two models in the folded Brillouin zone. For the five- 
orbital model, the results are obtained first in the un- 
folded Brillouin zone, and then those in the second Bril- 
louin zone of the original BCT lattice structure is moved 
to the first Brillouin zone. Here again they have simi- 
lar structures for most of the planes given here. The 
eigenvalue of the Eliashberg equation is 0.66 and 0.76 for 
s-wave, and 0.51 and 0.60 for c?-wave for the ten- and 
five-orbital models, respectively. Since the Fermi surface 
is different between the two models for some portions, 
the gap function differs between the two at those portions 
of the Fermi surface. Nonetheless, the present compari- 
son shows that the five-orbital model can be considered as 
reliable as far as the form of the superconducting gap on 
most of the k^ planes are concerned. Hereafter, we con- 
centrate on the five-orbital model and go down to lower 
temperatures, where larger-number fc-point meshes and 
Matsubara frequencies are required. 
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FIG. 4: (Color online) Comparison of the s± gap function 
between (a) the ten-orbital and (b) the five-orbital models at 
T = 0.07eV with 16 x 16 x 16 fc-point meshes. Contour plots 
are shown for five horizontal cuts around F-Z and X-P-X in 
the original Brillouin zone of the BCT lattice. The five-orbital 
results, obtained originally in the unfolded Brillouin zone, are 
refolded into this Brillouin zone. Solid black lines are the 
Fermi surface and the red dashed lines are the nodal lines of 
the gap. The Brillouin zone is shown along with the regions 
where the gap is presented. 
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(a)n=5.5 (b)n=5.8 (c)n=6.1 
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FIG. 5: (Color online) The spin susceptibility for the 
five-orbital models at T = 0.04eV with (a) n = 5.5 
(KFe2As2), (b) n = 5.8 (Bao.6Ko.4Fe2As2), and (c) n = 6.1 
(Ba(Feo.9Coo.i)2As2). (d) The Stoner factor against the band 
filling for the models of KFe2As2 and BaFe2As2. The dots 
indicate the positions of the actual band filling of the cor- 
responding material. At n = 5.8, we adopt the model of 
Bao.6Ko.4Fe2As2. 



V. FIVE-ORBITAL MODEL AT LOWER 
TEMPERATURE 

A. Spin fluctuation 

We now discuss the spin fluctuations using the five- 
orbital model. Here we take T = 0.04eV, adopting 
64 X 64 X 16 fc-point meshes and 512 Matsubara fre- 
quencies. One of the merits to adopting the five-orbital 
model is that it captures the essence of the spin fluc- 
tuations observed in neutron-scattering experiments. In 
FiglH we show the contour plot of the spin susceptibility 
of the five-orbital model in the unfolded Brillouin zone 
at kz — 0. In addition to the results for KFe2As2 (band 
filling n = 5.5), we also show for comparison the results 
for Bao.6Ko.4Fe2As2 {n = 5.8) and Ba(Feo.9Coo.i)2As2 
(n = 6.1). For Bao.gKo.4Fe2As2, we use the five-orbital 
model that is obtained by making a linear combination 
of the hopping integrals of BaFe2As2 and KFe2As2 us- 
ing the lattice structure of Bao.6Ko.4Fe2As2 given in Ref. 
ISlI as was done in Ref. l32l . For n — 6.1, we simply 
use the five-orbital model of BaFe2As2 adopting the ex- 
perimental lattice structure. '^^ The interaction reducing 
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FIG. 6: (Color online) The orbital character distribution on 
the Fermi surface for — 0. xy (a),(c) and yz (b),(d) or- 
bitals for the band fillings n = 5.5 (a),(b) and n = 6.1 (c),(d). 
The arrows indicate the wave vectors connecting the portions 
of the Fermi surface having similar orbital character. These 
arrows correspond to the peak position in the spin suscepti- 
bility. 



rate is taken to be / = 0.40, so that the Stoner factor 
does not exceed unity in all of the three cases considered 
here. In Fig. [SJa), the peak position is located at an in- 
commensurate wave vector [(l-2(5)7r,0], S = 0.17, which is 
in very good agreement with a recent neutron-scattering 
experiment for KFe2As2^. 

For Bao.6Ko.4Fe2As2 {n = 5.8), the peak moves to 
the commensurate position (7r,0). In the case of elec- 
tron doping, n = 6.1, the peak moves toward (tt, tt) 
in the unfolded Brillouin zone, which is consistent with 
the experimental results of Ba(Fe,Co)2As2 in the high- 
energy regime. The above tendency has also been ob- 
tained in other theoretical studies that use realistic band 
structures 

In the above, the adopted band structure is different 
among Figs. [5Ua),5(b) and 5(c). To focus on the evolu- 
tion of the spin fluctuations, we have also flxed the band 
structure to that of KFe2As2 or BaFe2As2, and varied 
the band flUing hypothetically from 5.5 to 6.2. We find 
that for both of the band structures, the spin fluctua- 
tion becomes commensurate around n — 5.8 ^ 6.0. For 
both of the models, it is found that the peak position 
continuously moves from the commensurate to the in- 
commensurate position as holes are dopedj^S, indicating 
that the spin fluctuation has the same origin in the entire 
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(Ba,K)Fe2As2 system, namely, the interband scattering 
between the band around (0,0) and (tt, 0)/(0, tt). 

The Stoner factor obtained for the above models is 
shown in Fig. EJd) as functions of the band filling. 
This shows that the strength of the spin fluctuations is 
stronger for the model of BaFe2As2 than in KFe2As2 
when compared at the same band filling. Therefore, 
when the composition of the material is changed continu- 
ously from Ba to K, we would expect the spin fluctuation 
to be maximum around the undoped regime (the Ba- 
rich regime), and decreases when the holes are doped by 
Ba— replacement. However, in the heavily hole doped 
regime (the K-rich regime) , the strength of the spin fluc- 
tuations should barely decrease, and may become even 
stronger with hole doping approaching the end material 
KFe2As2. 

As for the electron-hole asymmetry of the incommen- 
surability of the spin fluctuations, we have analyzed its 
origin in Ref. 23 as follows. The spin fluctuations are 
mainly governed by the interaction between portions of 
the Fermi surface having similar orbital character, so 
in Fig. IHa)-(d), we show the strength of the orbital 
character on the Fermi surfaces of KFe2As2 {n = 5.5) 
and BaFe2As2 with n = 6.1 in the unfolded Brillouin 
zone. Here we plot the states within a finite energy range 
-AE < E{k)-EF < AE with AE = 0.02eV, which can 
contribute to the spin fluctuations. Although Fermi sur- 
faces near the wave vector (tt, 0)/(0, tt) are barely present 
in KFe2As2, there are states in the vicinity of the Fermi 
level originating from a nearly flat band lying close to the 
Fermi level. The position of this flat band with respect to 
the Fermi level can be clearly seen in the magnifications 
of Figs. [TJa) and 1(b). As mentioned in Ref. [13, the 
spin fluctuations develop at wave vectors which bridge 
the portions of the Fermi surface having similar orbital 
character. The wave vectors connecting xy orbital por- 
tions and those connecting yz (or xz, not shown) portions 
nearly coincide in each of the two cases, so that this wave 
vector should correspond to the peak position of the spin 
fluctuations. In the electron-doped case, the connecting 
wave vector deviates from (7r,0) toward (tt, tt), while in 
the hole-doped case it deviates toward (0,0). Thus, the 
difference of the spin fluctuation incommensurability be- 
tween electron- and hole-doped cases originates from the 
multiorbital nature of the system. 



B. Superconducting gap 

Finally, we discuss the superconducting gap of 
KFe2As2. We obtain the superconducting gap at T = 
0.04eV for s-wave and dx2_y2--wave pairings. Here we 
take 32 X 32 X 16 fc-point meshes and 1024 Matsubara fre- 
quencies. We show in Fig. [7] the s-wave (a) and d-wave 
(b) gap functions of the flve-orbital model refolded into 
the original folded Brillouin zone. Here we take the in- 
teraction reducing rate / — 0.49. The rate / is increased 
compared to the values adopted for the spin suscepti- 



bility because here we focus only on KFe2As2, so that 
the Stoner factor does not exceed unity even for this /. 
We take larger values of the interaction so as to enhance 
the eigenvalue of the Eliashberg equation, which enables 
us to obtain a more reliable form of the superconduct- 
ing gap. The eigenvalue of the Eliashberg equation is 
As ^ 1.01 and Xd^2^y2 0.98, close to unity. These 
values indicate that the pairing symmetry competition 
is close. In fact, the possibility of d-wave pairing was 
proposed in a recent functional renormalization group 
study2i and an effective Hamiltonian approach. Even 
at this low temperature, the gap is essentially similar to 
the one obtained at higher temperature in Sec. IIVI The 
s-wave gap is basically s±-wave in the sense that the gap 
function changes sign between the states near k ~ (0, 0) 
and k ^ (tt, 0)/(0, tt). The origin of this sign change 
is strongly related to the spin fluctuation discussed in 
the present section. Namely, although the Fermi sur- 
face is barely present around k ^ (tt, 0)/(0, tt), there 
are states in the vicinity of the Fermi level originating 
from the bands that produces the electron Fermi surface 
for smaller hole content. The pair scattering from those 
states to the hole Fermi surface and vice versa mediated 
by the incommensurate spin fluctuations is the origin of 
the sign change of the gap function. On the other hand, 
the d-wave pairing originates essentially due to the in- 
traband repulsive interaction. In fact, we flnd that the 
spin fluctuations originating from the intraband interac- 
tion (within the hole Fermi surfaces) develops especially 
for the band filling below n — 5.5. 

Although the origin of the stabilization of the s-wave 
gap is the incommensurate spin fluctuation and the gap 
function (eigenfunction of the Eliashberg equation, to be 
strict) changes its sign between k ~ (0, 0) and k ^ (tt, 0), 
theoretically, the Fermi surface around the latter wave 
vector is barely present for KFe2As2, so, experimentally, 
this sign change should not be detectable. However, there 
are several possibilities of (nearly) vanishing gap func- 
tion, which may be relevant to the experimental obser- 
vations that suggest the presence of nodes in the super- 
conducting gap. First, the s-wave superconducting gap 
on the xy Fermi surface is rather small. One can see 
in Fig. [Tja) that the nodal line runs somewhat close 
to the xy Fermi surface. This may experimentally be 
detected as nearly vanishing gap, although there is no 
actual sign change in the gap function here. Second, the 
s-wave gap on the ai Fermi surface has different signs 
between the kz plane that includes the F point and that 
with the Z point, as seen in Fig. [Tl^a). This means that 
the gap has horizontal nodes on the ai Fermi surface, 
which is schematically shown in Fig. [Tjc). Although 
the five-orbital model is not reliable concerning this por- 
tion of the Fermi surface, this sign change is also found 
in the ten-orbital model (at higher temperature) as seen 
in Fig. |4] in Sec lIVI The sign change in this ai hole 
Fermi surface is a common feature found in the RPA 
studies of the 122 systems^^i^i and indeed some recent 
experiments2£ii£ support this possibility in phosphorous 
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FIG. 7: (Color online) The contour plots of gap function of 
the five-orbital model for s-wave (a) and cij.2_y2-wave (b) at 
T — 0.04eV. The solid black lines represent the Fermi sur- 
faces, while dashed red lines are the nodes of the gap. (c) 
A schematic figure of the horizontal nodes in the ai Fermi 
surface in the case of s-wave pairing. The arrow indicates the 
intraband scattering that drives the sign change. 



doped 122 materials. Third, there is a possibility of a 
d-wave gap, which is in close competition with the s- 
wave. This possibility, however, may not be consistent 
with an experiment on the magnetic vortex lattice, which 
shows that the gap is not strongly anisotropic within the 
planes.— 



VI. CONCLUSIONS 

In the present study, we have constructed an approxi- 
mate five-orbital model of KFe2As2 and compared it with 
the original ten-orbital model. We find the former to be 
reliable as far as the peak position of the spin fluctuations 
and form of the gap function are concerned. Using the 
five-orbital model, we have investigated the spin fluctua- 
tions and the spin-fluctuation-mediated superconducting 
gap of KFe2As2. We find that the spin fluctuation has an 
incommensurate peak structure, whose origin is the mul- 
tiorbital nature of the Fermi surfaces. The peak position 
of the spin fluctuation is in excellent agreement with the 
experiment 
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As for the superconductivity, s-wave and d-wave pair- 
ings are found to be in close competition. The s-wave 
gap has essentially the s±-wave form in the sense that 
the gap function changes its sign between the wave vec- 
tors around (0,0) and (tt, 0)/(0, tt). Although there are 
barely Fermi surfaces around (tt, 0)/(0, tt), the states in 
the vicinity of the Fermi level are nonetheless effective 
to give the incommensurate spin fluctuations and thus 
the sign change in the superconducting gap function. 
Although the sign change exists in the gap function, 
this should not be detectable in experiments because the 
Fermi surface is barely present around (0, 7r)/(7r, 0). As 
for the explanation for the nodes in the superconducting 
gap observed experimentally, besides a possible d-wave 
pairing, we raise two possibilities for the s-wave pairing: 
either horizontal nodes in the xzjyzj Fermi surface or 
the small gap on the entire xy Fermi surface. 
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